Water permeation through stratum corneum lipid bilayers from atomistic simulations 
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Stratum corneum, the outermost layer of skin, consists of keratin filled rigid non-viable corneocyte cells sur- 
rounded by multilayers of lipids. The lipid layer is responsible for the barrier properties of the skin. We calculate 
the excess chemical potential and diffusivity of water as a function of depth in lipid bilayers with compositions 
representative of the stratum corneum using atomistic molecular dynamics simulations. The maximum in the 
excess free energy of water inside the lipid bilayers is found to be twice that of water in phospholipid bilayers 
at the same temperature. Permeability, which decreases exponentially with the free energy barrier, is reduced 
by several orders of magnitude as compared to with phospholipid bilayers. The average time it takes for a water 
molecule to cross the bilayer is calculated by solving the Smoluchowski equation in presence of the free energy 
barrier. For a bilayer composed of a 2:2:1 molar ratio of ceramide NS 24:0, cholesterol and free fatty acid 24:0 
at 300 K, we estimate the permeability P = 3.7 x 10~ 9 cm/s and the average crossing time T m , = 0.69ms. The 
permeability is about 30 times smaller than existing experimental results on mammalian skin sections. 



PACS numbers: 82.70.Uv, 87.10.Tf, 87.15.Vv 

I. INTRODUCTION 

The 10 — 40 jiim thick stratum corneum (SC), the outermost 
layer of the skin |1], comprises rigid keratin filled pancake 
shaped non-viable cells (corneocytes) in an extra-cellular lipid 
matrix. This SC lipid matrix is believed to be the main barrier 
against water loss, with permeability values three to five or- 
ders of magnitude smaller compared to phospholipid bilayers 
forming plasma membranes Q]. This high permeation barrier 
is vital for life to maintain the required ionic concentration in- 
side the body for proper biochemical reactions. At the same 
time, controlled modulation of the barrier properties has the 
potential for non-invasive drug delivery Bfland restoration of 
the barrier function in compromised skin 14J]. 

The SC lipid matrix is conspicuous in having a large frac- 
tion of lipids from the ceramide (CER) family with long and 
asymmetric acyl tails. The other major components of the 
lipid matrix are cholesterol (CHOL) and free fatty acid (FFA) 
J5l0]. The bricks and mortar model 0], one of the widely ac- 
cepted models for the SC arrangement, pictures corneocytes 
as essentially impermeable bricks joined by the thin lipid ma- 
trix forming the mortar phase. One of the explanations put 
forward for the extremely low permeability of the SC is that 
the permeating molecules traverse a tortuous route through the 
lipid matrix. In this picture, the special chemical structure of 
the lipids plays no important role for the passive permeation. 
An apparent experimental justification for this picture is a very 
large lag time between introducing a radioactive molecule at 
one side of the SC and detecting it on the other side. With a 
simple diffusion model, even for water, one needs to invoke 
a path length two orders of magnitude larger than the phys- 
ical thickness of the SC layer considered @]. This picture 
assumes that water does not penetrate the corneocytes. How- 
ever, the corneocytes contain small hygroscopic molecules 



such as amino acids J^] (collectively referred to as natural 
moisturizing factor or NMF iflolfTlll ). For fully hydrated SC, 
water diffusivity in the corneocyte is estimated to be within a 
factor 2-3 of the diffusivity in bulk water lfl^[l3ll . So, it is a 
bit contradictory to consider them as acting as an impenetrable 
barrier against water transport. 

In this work, we limit our studies to simulations of hydrated 
lipid bilayers alone but include the minimal chemical details 
appropriate for the SC lipid matrix, and perform a series of 
molecular dynamics simulations to probe the permeability of 
the fully hydrated bilayers. Ceramide sphingolipids contain 
a fatty acid tail attached to a sphingosine motif. The fatty 
acid tail is highly polydisperse in length fTill . Also, there 
are at least 9 different classes of ceramides in human stra- 
tum corneum, with slight variations of the head groups, and, 
in case of ceramide 1, an additional esterified long fatty acid 
attached to the longer hydrocarbon tail. Free fatty acids also 
have a large polydispersity in vivo lfl5ll . For simplicity, we 
only consider ceramide NS 24:0 (ceramide 2), with its fatty 
acid tail containing 24 carbons. Similarly, the only free fatty 
acid we consider contains 24 carbon atoms. This particu- 
lar length was chosen to represent the majority fraction of 
ceramide and free fatty acid present in the human stratum 
corneum lfl4l[l5ll. 

CER with their long asymmetric tails control the main 
distinguishing features of SC lipid membranes compared to 
phospholipids lfl6ll . The absence of any large head group leads 
to close packing of hydrocarbon tails. We have performed 
simulations of pure CER bilayers, along with more realistic 
bilayers with CER, CHOL and FFA present in either 1:1:1 or 
2:2: 1 molar ratios. Simulation studies of water permeation for 
phospholipid bilayers, using very similar force-field as em- 
ployed in this study, exist in literature at 350 K lfl7tl . To com- 
pare and contrast between phospholipid bilayers and SC lipid 
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bilayer, we also simulated at 350K. For the 2:2:1 composition 
ratio, where comparisons are made with experimental results, 
we have done additional simulations at 300 K. 

Our main findings are that the diffusion of water molecules 
inside the lipid bilayer is highly anisotropic, and the ex- 
cess chemical potential for water is much higher than in 
typical fluid phospholipid bilayers at the same temperature. 
In the presence of the high free energy barrier, the simple 
absorption-diffusion picture, with which much of the exper- 
imental results are interpreted in the literature, is no longer 
valid. Instead, the time taken by a water molecule to cross 
the bilayer is determined by the Kramers' first passage time 
across the bilayer. Our results suggest that for water, the ex- 
perimentally determined lag time and thickness of the stratum 
corneum can be reconciled without invoking concepts like tor- 
tuosity when the free energy barrier is accounted for correctly. 



II. COMPUTATIONAL DETAILS 

To calculate the permeability coefficient from simulations, 
we constrain a water molecule at a fixed z distance from the 
lipid bilayer mid-plane. Here, the z direction is the direction 
normal to the bilayer with the bilayer center of mass being at 
z = 0. In this paper, we use the subscript _L to denote quanti- 
ties in the x — y plane of the bilayer. The average z component 
of the force on the constrained molecule is related to the spa- 
tial derivative of free energy dill 



dAG(z) 
dz 



(1) 



Here, (• • • ) refer to averages both over time and different sys- 
tem replicas. The value of the local excess chemical potential 
can be calculated through numerical integration from the bulk 
water phase. The auto-correlation of the force is related to the 
diffusion coefficient through a Kubo relation 0T 
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where, R is the universal gas constant, T the absolute tem- 
perature and AF z (z,t) = F z (z,t) — (F z (z)). The macroscopic 
permeability coefficient P, the ratio between the flux and con- 
centration difference (J = PAc), is macroscopically defined 
through HI 
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where the integration runs over the bilayer thickness, = 
ksT, with kg being the Boltzmann's constant. The constraint 
on the water molecule is only along the z direction, so the in- 
plane diffusivity D± in the x — y plane is calculated from the 
mean square displacement as a function of time, (Arj_) (?) = 
4D ± t. 

To estimate x av , the average time it takes for a water 
molecule to cross a bilayer, we introduce a steady current J 
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FIG. 1: Consideration for calculation of the average time l av for 
crossing the bilayer. At one side of the bilayer a steady current is 
supplied, which is forced to go through the bilayer because of a hard 
wall on the other side of the insertion point. At the maximum of the 
excess chemical potential an absorbing wall removes any molecule 
reaching that point. 



at z = —d, at one side of the bilayer (Fig.[TJ. The number den- 
sity of water molecules n(z) obeys the Smoluchowski equa- 
tion (H, 

dn d n (dn . dAG\ TS . J , 



We consider a hard reflecting wall at z < —d, where = o 
(in the bulk AG = 0). An absorbing wall is placed at the max- 
imum of the potential z = Zm- Once steady state is achieved, 
the total number of surviving molecules is 



n(z)dz = J% a 



(5) 



In the steady state =0), integrating Eq.|H for z > d leads 
to 



n{z) = Je 
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Using this expression in Eq.[5]gives 
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We compute x av by numerical integration of functional fits to 
AG(z) andD(z). 

All simulations are done with extended ensemble molecular 
dynamics at constant temperature and pressure ensemble with 
GROMACS molecular dynamics package JM, HH IH with a 
timestep of lfs. The interaction parameters are based on the 
united atom OPLS force field Q23h with modifications for the 
nonpolar hydrocarbon groups 12411 that accurately reproduce 
experimental results for lipid molecules l25ll . Polar hydrogens 
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were included explicitly. The dihedral potentials in the lipid 
tails were described by the Ryckaert-Bellemans term j26ll . At 
skin conditions the SC lipids do not fully ionize, so, polar 
groups were assigned partial charges chosen from previous 
simulations of similar molecules 11271 l28ll . The SPC model 
ll29ll was used to describe the water molecules. 

Nose-Hoover thermostats ll30ll3Tll coupled separately to the 
lipid and water molecules with a time constant of 5ps were 
used to control the temperature. The Parrinello-Rahmanbaro- 
stat ll32l l33ll with time constant 5ps was used for pressure 
coupling. The diagonal components of the compressibility 
matrix were chosen to be 4.5 x 10~ 5 bar. The off-diagonal 
components were set to zero to keep the simulation box or- 
thogonal and standard periodic boundary conditions were ap- 
plied in all three directions. Electrostatic interactions were 
calculated with a group-based cut-off. The cut-offs for both 
the Van der Waals and electrostatic interactions were set to 
1.2nm. With the small dipole moments involved in these 
simulations, the electrostatic interaction becomes negligible at 
sufficiently small distance so that the results remain indepen- 
dent of using either a group-based cut-off, or Ewald summa- 
tion to take account of interactions with the periodic images 
llrjll . All lipid bonds were constrained with the SHAKE algo- 
rithm l34ll . Rigid SPC water molecules were updated with the 
analytic SETTLE algorithm JH. 

From the final equilibrated configurations from a previ- 
ous study UH], we selected a random water molecule ap- 
proximately 5nm above the bilayer midplane. The water 
molecule was pulled along the negative z direction at a rate 
of 0.05nm/ps by moving the water molecule 5 x 10~ 5 nm ev- 
ery (discrete) fs time step. Each time the relative z separa- 
tion between the water molecule and the bilayer center of 
mass changed by 0.2nm, the configuration was equilibrated 
for 200 ps with the z separation between the water molecule 
and the bilayer midplane kept fixed. The configurations at 
these steps were stored at full precision for further calcula- 
tions. 

We evolved each of these saved configurations, with the se- 
lected water molecule constrained to be at fixed z separation 
from the bilayer midplane for 2 ns. At each time step, the force 
F z along the z-direction on the constrained water molecule and 
the in-plane displacement of the center of mass of the con- 
strained water molecule were stored. Every lOOps we also 
store the indices of the atoms within a distance of 0.4nm to 
find out about the local environment of the constrained water 
molecule. The whole procedure was repeated with 15 differ- 
ent random water molecules for each composition and tem- 
perature investigated. 

To calculate D z from the autocorrelation of F z (Eq. [2]), we 
first estimate the decay time T e of < F z (t)F z (0) >. Typical val- 
ues of T (: are ~ 0.1 ps. The upper time limit of the integration 
(eqn.|2]i to calculate D z was chosen to be 100t c . We calculate 
the in-plane diffusivity D± from long time behavior (> 0.5 ns 
) of the mean-squared displacements with time origins chosen 
at intervals of 1 ps for averaging. 

Data from simulations with different water molecules were 
used together to calculate the average and error estimates for 
F z (z), D z (z) and D±(z). Because the bilayers studied here 



are symmetric, D z (F z ) is an even (odd) function of z- Using 
the data from both leaflets we achieve a better estimate at a 
given z- Numerical integration of F z (z) from the bulk gives 
the excess free energy AG(z) (Eq.Q~|i. Finally, permeability is 
calculated by numerical integration across the bilayer (Eq.O 
using the calculated D z (z) and AG(z). 



III. PERMEABILITY OF CERAMIDE BILAYER 

The asymmetric long chain ceramides are responsible for 
much of the distinguishing features of SC lipid structure as 
compared to other biologically relevant membranes. In this 
first part of this work we concentrate on a pure CER bilayer 
and contrast our findings with literature results on phospho- 
lipid bilayers. 




FIG. 2: Snapshot of hydrated CER bilayer at 350 K with a con- 
strained water molecule shown as a large sphere. 

Fig. [2] shows a snapshot of CER bilayer containing 128 
CER molecules and 5250 water molecules with a water 
molecules constrained to be at a distance 1.05nm from the 
bilayer midplane along the z direction. The long hydrocarbon 
tails lead to large nematic order in the bilayer. Once the wa- 
ter molecule is inside the lipid layer, it faces little resistance 
in moving along the z direction, but the motion in the x — y 
plane is severely restricted. Because the two tails have large 
asymmetry, the midplane region is mostly occupied by atoms 
from the longer tail of ceramide and there the nematic order is 
lower. 

The diffusivity of water as a function of distance from 
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FIG. 3: Diffusion coefficients of water as a function of distance from 
the bilayer midplane for hydrated CER bilayer at 350 K. The diffu- 
sivity along the normal direction to the bilayer, D z , is shown as filled 
square and the in-plane diffusivity in the x — y plane, Dj_ is shown as 
empty circles. The position of the head group and the two tails are 
indicated by the schematic representation of the ceramide molecule 
below the plot. We exploit the symmetry about the bilayer mid-plane 
to get better averages. Thus data on only one side of z = is inde- 
pendent. 



the bilayer midplane is shown in Fig. [3] In the bulk water 
(|z| > 4nm), both methods lead to very similar value of the 
diffusivity. D± falls by nearly three orders of magnitude in the 
ordered tail region, while D z is reduced by slightly less than 
an order. This agrees with the interpretation of water molecule 
inside the ordered region of the bilayer being effectively con- 
strained in channels defined by the hydrocarbon tails. Close to 
|z| ~ 2nm, the measured D± from the simulations fall below 
10~ 6 cm 2 /s. There, the water remains confined in the same 
lipid neighborhood for the entire simulation and we can no 
longer calculate D± with certainty from the mean square dis- 
placement. Except for this narrow region (\z\ — 2nm), the 
diffusivity is high enough for the water molecule to explore a 
large part of the system in the x — y direction within the run 
time of a single simulation. The asymmetry between the two 
tails lead to a low density molten region at the bilayer mid- 
plane, where both D z and D± approach the bulk diffusivity 
and is isotropic. The decrease in the value of D z inside the 
bilayer is comparable to that found in DPPC bilayer simula- 
tions B17I1 . However, in DPPC, D± remains comparable to 
D z throughout the bilayer. At the midplane of DPPC bilayer, 
because of large free volume, the diffusivity was found to be 
almost twice that of bulk water. In ceramide, both D z and D± 
approach the bulk water value at the midplane, but remain less 
than it. The asymmetry in the two tail lengths allow the CER 



bilayer to have a much larger local density at the midplane 
region than the DPPC bilayer, because of partial interdigita- 
tion. This limits the diffusivity to a value lower than that in 
the bulk. 
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FIG. 4: Excess chemical potential for a water molecule inside a CER 
bilayer at a given depth at 350 K. z = is the bilayer mid-plane. 

From the average force measured as a function of the depth, 
we numerically compute the excess chemical potential AG of 
a water molecule as a function of z (Fig. |4j. The chemical 
potential rises steeply on entering the lipid bilayer, reaching a 
maximum (~ 43kJ/mol ~ l5ksT) at a distance of ~ 1.5nm 
from the bilayer midplane. At the midplane itself, the excess 
chemical potential drops to a value of ~ 29kJ/mol(~ lOksT). 
The maximum in the excess chemical potential is about twice 
as high as that measured from DPPC bilayer simulation at the 
same temperature |fl7ll . 

The permeability decreases exponentially with the excess 
chemical potential. Our estimate for the permeability for CER 
bilayer at 350K, P = 1.12(±0.06) x l(T 8 cm/s, is approxi- 
mately five orders of magnitude smaller than the DPPC bi- 
layer permeability of 7(±3) x 10~ 2 cm/s lH7ll at the same tem- 
perature. This is mainly due to the much larger free energy 
barrier and partially due to reduced diffusivity compared to 
DPPC bilayer. 

Using the excess free energy and diffusivity, numerical in- 
tegration of Eq. [7] leads to the average crossing time for ce- 
ramide bilayer x av = 3.3 x 10~ 4 s. 



IV. MODEL STRATUM CORNEUM LIPID BILAYERS 

To provide some understanding for experiments with bi- 
ological skin sections, we consider two ternary mixtures of 
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FIG. 5: Snapshots of bilayers at 350K containing a 2:2:1 molar ratio 
of CER, CHOL and FFA. Part of the lipids were stripped along the 
x axis to show the molecular arrangement near the water molecules 
(shown as large spheres). The CER, CHOL and FFA molecules are 
colored as cyan, green and maroon respectively. In the top panel, 
only the constrained water enters the bilayer and in this particu- 
lar frame is contained in the free volume created by shorter CHOL 
molecules. In the bottom panel the constrained water molecule facil- 
itates entry of another water molecule in the lipid bilayer. 



CER, CHOL and FFA with 1:1:1 and 2:2:1 molar ratios. The 
2:2:1 composition is considered to be representative of the 
in vivo composition y]]. For this composition, besides 350 K, 
we perform simulations at 300 K, close to the physiological 
temperature. In this section we concentrate on the 2:2:1 mix- 
ture and consider the results on the 1:1:1 mixture in Table U 

CHOL reduces the local nematic order in the three compo- 
nent bilayers. Also, to fit the rather bulky CHOL molecule, the 
bilayer needs to create more free volume. This effect on the 
free volume is opposite to that in phospholipid, where choles- 
terol was found to reduce the available free volume in DPPC 
bilayer I3r3l . In the simulations with CER bilayers only the 
constrained water molecule enters the bilayer, while in the 
three component system another water molecule occasionally 
joins the constrained water molecule (Fig. |5). This reduces 
the energy due to favorable hydrogen bonds between the two 
water molecules. At the same time, especially close to the 
CHOL molecules, large enough free volume is available to 
accommodate the two water molecules. 
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FIG. 6: Top panel: Diffusion coefficients D z (open symbols) and Dj_ 
(filled symbols) for water molecule as a function of distance from the 
bilayer midplane of containing a 2:2:1 molar ratio of CER, CHOL 
and FFA respectively at 300 K (circles) and at 350 K (squares). Bot- 
tom panel: Mass densities as a function of z at 300 K. 

Fig. [6] shows the water diffusivity as a function of z for the 
2:2: 1 bilayer. At the bottom of the figure, we show the local 
mass density of the three components at 300 K. FFA follows 
the CER distribution closely, while CHOL prefers to stay just 
below the CER head group, finding favorable hydrogen bond- 
ing with CER. Water diffusivity in the bulk scales normally 
with temperature, varying from 7.5 x 10~ 5 cm 2 /s at 350K to 
4.5 x 10~ 5 cm 2 /s at 300 K fH. At 300 K, D z and D ± approach 
each other only near \z\ ~ 5nm, signifying considerable or- 
dering of the water molecules in the bulk liquid close to the 
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T 


molar ratio 




AGmax 


pAG max 


P 






(K) 


CER:CHOL:FFA 


(10~ 5 cm 2 /s) 


(kJ/mol) 




(cm/s) 




(ms) 


350 


1:0:0 


1.3 


42.7 


14.7 


1.1 x 10" 


s 


0.33 


350 


1:1:1 


1.1 


34.6 


11.8 


8.2 x 10" 


-8 


0.03 


350 


2:2:1 


1.0 


33.8 


11.6 


1.3 x 10" 


-7 


0.02 


300 


2:2:1 


1.1 


38.5 


15.3 


3.7 x 10" 


-9 


0.69 



TABLE I: Summary of the main results from the simulations. 



bilayer. Inside the bilayer the diffusivity is much less affected 
by the temperature, because the tail ordering and free volume 
do not change by much for this system in this temperature 
range. The range of z, in which the water molecule is es- 
sentially confined in the same lipid neighborhood during the 
entire simulation, is larger than for the pure CER bilayer. 
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FIG. 7: Excess chemical potential for a water molecule inside a bi- 
layer containing a 2:2:1 molar ratio of CER, CHOL and FFA. 

The maximum in the excess chemical potential (Fig. [7]l is 
lower than in the pure CER bilayer (34kJ/mol at 350 K and 
38kJ/mol at 300 K). At 350K, the permeability of 2:2:1 bi- 
layer P = 1.30(±0.01) x 10" 7 cm/s, is about one order higher 
compared to CER bilayer. Permeability drops by nearly two 
orders of magnitude on reducing the temperature to 300 K, 
giving, P = 3.67(±0.02) x 10~ 9 cm/s. The average crossing 
time for the 2:2: 1 bilayer increases from 0.2 x 10~ 4 sat350K 
to 6.9 x 10~ 4 s at 300 K. 



2:2:1 bilayer at 350K is ~ 0.7g/cc IToTl . comparable to that of 
liquid hexadecane ~ 0.73 g/cc [13811 . The excess free energy of 
water in hexadecane is ~ 25kJ/mol l39Tl . which agrees with 
the excess chemical potential at the bilayer midplane (Fig. |7j. 
The strong nematic order induces much larger density close 
to the head groups, with average density for the SC bilayers 
being ~ 0.94 g/cc 111611 . This reflects in the maximum in free 
energy being ~ 40kJ/mol. Arrhenius plots of the temperature- 
dependent permeability from human [2] and porcine 14011 stra- 
tum corneum suggest an activation energy 60kJ/mol. 

Experiments on human ll4ltl and porcine [ 8] SC found a per- 
meability of order 10~ 7 cm/s. The permeability value at 300K 
for the 2:2:1 bilayer from our simulations is 3.7 x 10~ 9 cm/s, 
which is about a factor of 30 lower than the experimental val- 
ues. With exponential dependence of AG, permeability de- 
pends strongly on the composition considered. In our simula- 
tions, we found that the 2:2:1 bilayer shows one order higher 
permeability compared to the CER bilayer at the same tem- 
perature. The lipids in the stratum corneum are highly poly- 
disperse and contain unsaturated lipid tails. The presence of 
unsaturated fatty acids and polydispersity will probably in- 
troduce greater disorder, reducing the permeability as com- 
pared to the 2:2:1 bilayer considered in this study. Also, a 
patch of stratum corneum is unlikely to have a defect free 
lipid structure throughout the sample. Confocal laser scan- 
ning microscopy seems to suggests that the stratum corneum 
shows large variability in permeation over fim length scales 
1I42L If there are defect pathways that offer less resistance 
than the defect free lipid layers, the experimental result on the 
permeability is likely to be dominated by these defects. Both 
the increased disorder due to molecular polydispersity and the 
presence of defects in the macroscopic sample will increase 
the permeability above that of a single perfect bilayer (as sim- 
ulated). 

In permeation experiments there is a delay between the in- 
troduction of radioactive water vapour at one side of the sam- 
ple and its first detection on the other side. This lag-time T/ aA , 
was experimentally found to be ~ 8 \ 10 3 s for porcine skin at 
303 K|@]. In a Fickian diffusion model, this would correspond 
to a diffusion path length 8 ~ ^jD%\ af ,. Reasonable values of 
the diffusivity D lead to a path length which is much larger 
than the physical thickness of the sample JH], and hence the 
interpretation of a tortuous path that avoids the corneocytes. 

Concentrating on the 2:2:1 bilayer at 300 K, and using 
the bilayer thickness (~ 5nm) and the minimum in D z (~ 
10~ 5 cm 2 /s), the Fickian diffusion picture yields X[ ag ~ 10~ 8 s. 
This estimate is a factor of 7 x 10 4 lower than the calculated 
mean crossing time (6.9 x 10~ 4 s) which takes into account the 
free energy barrier (table UJ). Use of D z and X av in this Fick- 
ian diffusion picture will suggest an apparent pathlength that 
is ~ 250 times larger than the real bilayer thickness, although 
the molecule actually only traversed the bilayer thickness. 



V. COMPARISON WITH EXPERIMENTS 



VI. CONCLUSIONS 



In table U we summarize the main results for the CER bi- 
layer and the ternary mixtures. The midplane density for the 



In conclusion, we have calculated the excess chemical po- 
tential profile and diffusivity for water molecules in fully hy- 
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drated lipid bilayers composed of lipids corresponding to the 
stratum corneum lipid matrix. We find that, compared to phos- 
pholipids simulated with very similar force fields and at the 
same temperature as in this study, the SC bilayers show nearly 
twice as high a chemical potential barrier against water per- 
meation. The high degree of correlation of the free energy 
profile with the lipid density profile suggests that the main rea- 
son for this large free energy barrier is because the ceramide 
lipids lead to dense packing of the acyl tails. Water diffusiv- 
ity across the bilayer does not change drastically as compared 
to phospholipids. High degree of tail ordering ensures that 
once the water molecule is inside the hydrocarbon region, it 
can move without large extra free energy cost along the chain- 
orientation direction. However, diffusivity in the perpendicu- 
lar direction (in-plane diffusivity) is decreased by two orders 
of magnitude. 

The permeability coefficients from our simulations are 
much smaller than the experimental results. In our analy- 
sis, we showed that the high free energy barrier necessitates 
the use of a first passage time to estimate the lag time, and 
neglect of the barrier leads to an apparent pathlength that is 
much larger than the physical thickness of the sample. Our 



results suggest that, for water permeation, we do not need to 
invoke impermeable corneocytes. 

It is possible to discern the importance of the free energy 
barrier on the lag time from experiments at different temper- 
atures. The Kramers' mean first passage formalism suggests 
that the temperature dependence of the lag time will be Ar- 
rhenius like: ln(T/ a? ) ~ 1 fT, while a simple diffusive picture 
suggests Ti ag ~ l/T, We are not aware of any measurements 
which looked at this temperature dependence. 
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